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Abstract 

The primordial non-Gaussianity parameters /nl and tnl may be scale-dependent. We investigate the 
capability of future measurements of the CMB /i-distortion, which is very sensitive to small scales, and of 
the large-scale halo bias to test the running of local non-Gaussianity. We show that, for an experiment 
such as PIXIE, a measurement of the /i-temperature correlation can pin down the spectral indices n/ NL and 
n TNL to values of the order of 0.3 if /nl = 20 and tnl = 5000. A similar value can be achieved with an 
all-sky survey extending to redshift z ~ 1. In the particular case in which the two spectral indices are equal, 
as predicted in models where the cosmological perturbations are generated by a single field other than the 
infiaton, then the 1-cr error on the scale-dependence of the non-linearity parameters goes down to 0.2. 
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1 Introduction 



Detecting a possible primordial source of non-Gaussianity (NG) in the cosmological perturbations is one of the 
main targets of current and future experiments measuring the properties of the Cosmic Microwave Background 
(CMB) anisotropics and the large-scale structure. Indeed, measuring a certain level of NG in the three- (bispec- 
trum) and four-point (trispectrum) correlator of the perturbations opens up a unique window into the physics 
of inflation which is believed to be the period during which such fluctuations are quantum-mcchanically gener- 
ated [1] . The current constraints on NG come from the measurement of the CMB anisotropy bispectrum [2] and 
from the properties of the clustering of galaxies which has been identified to be a powerful probe of NG thanks 
to the fact that NG introduces a scale-dependent bias between the power spectra of halos and dark matter [3,4]. 

Most of the attention in the literature has been devoted to the so-called "local" model of NG, where the NG 
is defined in terms of the primordial gravitational potential $(x) as 



<f(f) = <j> G (x) + /nl [^(z) - (4(f)) 
The corresponding bispectrum and trispectrum are given by 

B*{k x ,k 2 M) = 2/ NL [> (fc 1 )iyfc 2 ) + 2 eye. 
25 r 

T^(ki,k 2 ,k 3 ,k 4: ) = yfNL P^{ki)P <t ,{k 2 )P < p{ki Z ) + 11 eye. 



(1) 

(2) 
(3) 



where P<j,(k) is the power spectrum of the gravitational potential. This type of NG is generated in multifield 
inflationary models where the cosmological perturbation is sourced by light scalar fields other than the inflaton. 
The corresponding perturbations are both scale invariant and special conformally invariant [5,6]. The parameter 
/nl is currently constrained to be in the range (32 ± 21) by WMAP [2] and (28 ± 23) by the large-scale 
structure [7], while the parameter r NL needs to be in the range (—0.6 < r NL /10 4 < 3.3) as inferred from the 
WMAP 5-year data [8] . Measuring the amplitudes of both the bispectrum and the trispectrum is extremely 
interesting as, if only one degree of freedom is responsible for the perturbations, then there is a well-defined 
relation between the NG parameters, tnl = (§/nl) • On the contrary, if more than one field is responsible 
for the cosmological perturbations generated through the inflationary dynamics, then there exists an inequality, 
t nl > (|/nl) [6,9,10]. To which extent future measurements of the scale-dependence of halo bias can test 
multi-field inequality has been the subject of Ref. [11]. 

Even though the definitions (2) and (3) are widely used to model NG in the primordial perturbations, it is 
just the first step one can make on this matter. One, more general, definition of the bispectrum and trispectrum 
could include a scale-dependence in the non-linearity parameters /nl and tnl- This step is well- motivated by 
the theoretical predictions of some models [12-16]. The running with physical scale of the NG parameters /nl 
and tnl has been the subject of an intense recent research [17-24]. 

To account for the running of /nl m its full generality one can adopt for example the parametrization used 
in Ref. [25] (see also Ref. [14]) 



B*(h,k 2> k 3 ) =2K An ,(fc 3 )f m (Ai)^ )n (fca)P^(fei)P^(*2) + eye. 



where 



(M 



fcr 



(4) 



(5) 



Here £/ NL (&) parametrizes the (self-)interactions of the fields and £ m (fc) the ratio of the contribution of each 
field. From this general parametrization, we can also easily extend the one for the trispectrum 



25 

T$(ki,k 2 ,k 3 , fc 4 ) = — [CT NTj (h,h)^m(h)^ m (k2)^m(h3)P4(ki)P,j,(k 2 )P^(ki3) + eye] , 



(6) 



where 



£t nl (ki, kj) — £ TNL (fco) 



kikj \ 



k 2 



(7) 



In the single-field limit, £, TNL (ki, kj) = §§£/ NL (&«)£/!«,(%) an d £m(k) = 1. According to this parametrization, in 
the case of a multi-field inflation, we have three free parameters, n/ NL , n m and fV NL , which describe the scale 
dependence of the non-linearity parameters /nl an d tnl and of the dimensionless power spectra. In order to 
decrease the complexity of the analysis, from now on we make the assumption that n rn is significantly much 
smaller than unity. By doing so, we are left with the following parametrization of the non-linear parameters 
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/ U \ ™/nl 

and 

TNL(ki,k j ) = T^( k ^\ . (9) 

CMB information alone, in the event of a significant detection of the NG component, corresponding to /nl = 50 
for the local model, is able to determine n/ NL with a Tcr uncertainty of about 0.1 for the Planck mission [17]. 
A local bias analysis performed in the same Ref. [17] showed that high-redshift surveys (z > 1) covering a 
large fraction of the sky corresponding to a volume of about 100 h~ 3 Gpc 3 might provide a 1-a error on the 
running /nl parameter of the order of 0.4(50//nl)- On the other hand, using the WMAP temperature maps, 
a constraint on the running of the scale-dependent parameter /nl has been recently obtained in Ref. [26] to be 
ii/ Nl = 0.30(+1.9)(— 1.2) at 95% confidence, marginalized over the amplitude / NL . To the best of our knowledge, 
no forecasts for the running of the trispectrum parameter tnl exist in the literature. In fact, in the case in 
which the perturbations are sourced by a single field, then a well-defined relation between the running spectral 
indices holds, 



"/nl = "r NL (10) 

and the indices are therefore not independent. In this paper we will assume that /nl and tnl, and therefore 
their spectral indices too, are not related to each other, thus leaving open the possibility that the perturbations 
are originated from a multi-field scenario. 

The goal of this paper is to provide some useful forecasts on the spectral indices n/ NL and ^ TNL from the 
possible physical imprints that NG can leave on the the CMB /i-distortion and the halo bias. Measurements 
of the /i-type distortion of the CMB spectrum provide the unique opportunity to probe these scales over the 
unexplored range from 50 to 10 Mpc~ 4 and it has been recently pointed out that correlations between fi- 
distortion and temperature anisotropics can be used to test Gaussianity at these very small scales. In particular 
the /x-temperature cross correlation is proportional to the very squeezed limit of the local primordial bispectrum 
and hence measures /nli while the fj,-fi is proportional to the primordial trispectrum and measures tnl [27] 
(see also [28]). Being the ^-distortion localized at small scales, we expect it to be very sensitive to the possible 
running of the NG parameters /nl and tnl- This will be studied in section 2. In section 3 we will study 
the effect of running NG parameters onto the halo bias, taking into account the running of the trispectrum 
amplitude as well. Our conclusions will be presented in section 4. In all illustrations, the cosmology is a flat 
ACDM Universe with normalisation er 8 = 0.803, Hubble constant h = 0.701 and matter content f2 m = 0.279. 



2 CMB /i-distortion 

The goal of this section is to compute the effect of the running NG onto the CMB /z-distortion. The latter is 
caused by the energy injection originated by the dissipation of acoustic waves through the Silk damping as they 
re-enter the horizon and start oscillating. The interesting property is that this effect is related to primordial 
perturbation scales of 50 < fcMpc < 10 4 which are not accessible from CMB anisotropics observations. 

At early times (z 3> z^j = 2 x 10 6 ), the content of the universe can be described by a photon-baryon fluid in 
thermal equilibrium which has a black-body spectrum. This equilibrium is achieved mainly through elastic and 
double Compton scattering. However, at later times (z^j = 5 x 10 4 < z < z^i), double Compton scattering 
is no longer efficient whereas the single Compton scattering still provides equilibrium. The photon number 
density is however frozen and only the frequency of the photons can be changed. It can be shown that any 
energy injection in the photon-baryon fluid will distort the spectrum by the creation of a chemical potential 
/i. The photon number density per frequency interval is then n{y) = (e xJrli ( x > — 1) _1 , where x = hv j '(ksT). 
The parameter fj, due to damping of acoustic waves can then be expressed in terms of the primordial power 
spectrum [29]. Using the Bose-Einstein distribution plus the fact that the total number of photons is constant, 
for an amount of energy (density) released into the plasma 8E/E, one finds that \i ~ 1A5E/E, where 

— ~- (5 2 Jx))\ z ^ , (11) 

and (5^(x)) represents the photon energy density fluctuation averaged over one period of the acoustic oscillations. 
As the modes of interest re-enter the horizon during the radiation phase, one finally finds 
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M (f) ~ 4.6 / ^/ C^e^l^ ) (co S (fc 1 r)cos(fc 2 r)) p e'^+^V^, ^ . (i 2 ) 



1 •-<■ 



where £(if) = 5 < f>(x)/3 describes curvature perturbations, fc± = k\ ± and in order to account for the fact that 
the distortion arises from a thermalization process, one uses a top-hat filter in real space W(x), which smears 
the dissipated energy over a volume of radius kp 1 ^ < k^ 1 , where k£>{z) is the diffusion momentum scale 

k D {z) ~ 4.1 • 10~ 6 (1 + z) 3/2 Mpc _1 . (13) 

We proceed by computing the correlations between /i-distortion and temperature anisotropy as well as fj.fi self 
correlation as done in [27], but allowing for a running of /nl and tnl given by Eq. (8) and (9). The curvature 
perturbation bispectrum in the squeezed limit (k 3 <C k\ ~ k2) is expressed as 

(C&C&Cjg,) = (27r) 3 5 3 (fci + fc 2 + fc 3 )^/ NL (fc_/2)P(fe_/2)P(fe+). (14) 
The temperature-// correlation therefore reads 1 

- -2.2X10--/^^, (15) 

where the primordial curvature spectrum is defined by (Cg C^ 2 ) = (27r) 3 (5 3 (fci + k)2ir 2 A 2 (fci) / fc 3 with A 2 (fcp) 
2.4 x 10 -9 at the pivot scale fc p = 0.002 Mpc -1 [2]. The parameter 6 is defined by 

f dlnk +n (k +n fw fk 



dlnfe_ 



V /NL A2(fc_/2)A2(fc + ) 



,-k 2 _/(2k 2 D (z)) 



(16) 



The /i-distortion is created during the period between z„ j = 2 x 10 6 and z„ * = 5x 10 4 which implies kp^. 
11600 Mpc -1 and k D j ~ 46 Mpc -1 . For a weak scale dependence A 2 (fc) = A^(fc p )(fc/A:p)™ s ~ 1 we obtain 



l i / i V s ' /' 1 

In f*R±\ n s + n /NL - 1 ^ y/2k p ) \ y/2K 



[knizr^^- 1 ]^. (17) 



If we take the same pivot for /nl as for the primordial spectrum, fc* = k p , the equation above becomes the 
same expression as for a constant /nl = /nl but with a shifted spectral index n s replaced by (n s + n / NL ) . This 
shows explicitly that we recover the scale invariant result for n/ NL = and we have b ~ 1 + 10(n s + n/ NL — 1) 
for (n a + n/ NL - 1) ~ 0. 

Using the trispectrum in the collapsed limit (k\2 ~ 0) 

(C&CfeCj^) = (2^) 3 <5 3 (fcx + k 2 + k 3 + fc 4 )4T NL (/c_/2, k 3 )P(k./2)P(k+)P(k 3 ), (18) 
we obtain the NG contribution to the //-distortion self-correlation 

- 5.6x10-^^, (19) 



lr To compute the temperature anisotropies we adopt the same approximation as in Ref. [27], that is the Sachs- Wolfe approxima- 
tion. Based on the findings in Ref. [28] , where the full radiation transfer function was adopted, we expect an overall decrease of the 
signal-to-noise ratio of order of 40%. Later in the text, we also point out that the change of the pivot scale amounts to corrections 
of the order of 30%. 
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where 



21(1 + 1) 



k-k 3 



din fc + dlnfc_ din fc 3 ji(k + r e ) 2 W 



k+ 



2kl 

1 



A*(k-/2)A*(k + )A}(k 3 ) 



2 ; 



' e -k 2 _/(2k 2 D )~ 




' e -2k 3 /(2k 2 D )~ 













■Tl.-l 



2(n s -l) 



2n TK 



[k D (zY 



(20) 



This is just & 2 with the index n / NL replaced by ?J rNL and it corresponds to the result of a constant tnl = t"nl 
with n s replaced by (n s + n rNL ). We recover the scale-invariant result for n TNL = 0. The behaviour of the 
parameters b and b is shown on Fig. 1. 
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Figure 1: Value of the parameters 6(n/ NL ) and b(n TNIj ) for two different pivot scales fc* and n s = 0.96. The dashed line 
shows the maximal value of n/ NL (and n rNL ), for which the approximation (25) is correct. 

Having computed the key parameters b and 6, we proceed by estimating the signal-to-noise ratio to estimate 
the values of n/ NL and jv nl measurable from the /^-distortion assuming that the amplitude and are 
known from some other experiments. In general the signal-to-noise ratio for variables Aj is defined in terms of 
the Fisher matrix as [30] 

| ee X \,F, J X ! . (21) 

In the case of only one variable, it reduces to S/N = \\/F = \/a\. In our case, to measure the spectral index 
"/nl we can adopt the Fisher matrix 




while for the spectral index n TNL we adopt the Fisher matrix 




The noise for /i-distortion can be modelled assuming a Gaussian beam experiment [31] by 

c^ N ~ w ;V 3 /^«, (24) 
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where £ max is the maximum multipole fixed by the experiment's beam size and is the sensitivity to /x. For 
the PIXIE experiment [32], £ max — 84 and = y/4ir x 10~ 8 . We also approximate the variance of the Ci's 

by 



2 



2f 



^ ((Cf + C^ N ){Cj T + Cj T < N ) + (Cff 



1 



2Z+1 



r TT r p.p, 



,iV 



and 



2 ,„„.„N\1 



21 + r 

where we used that Cj T » Cj T ' N , C^ N » and Cj T C^ N » (Cf^) 2 . 2 
The signal-to-noise for n/ NL at fixed is given by 

S \ 



A 



»/NL 



'«V In 



In 



/ 1 



ll^A^(fc p )/^ L 
1 



/ 1 



\V2h 



1 



n.—l 



V2~h p 



1 



n /NL + n 

whereas the signal-to-noise for n TNL at fixed r^ L is 



1 



n fNh +n s -l 
K D 



Mk D )k D f ^ 



+n s -l" 



307tw m A^(A;p)tnl In 2 



x In 



ln(fc£,)fc^ TI 



1," T NL 

D 



+n.-l 



6(n T 



(25) 
(26) 



(27) 



(28) 



The left plot of Fig. 2 shows ^/ nl (/nl) at (S/N)n f L = !• An amplitude < 10 2 enables to detect 
nj Nl > 0.3, at least with the PIXIE experiment. Notice also that the dependence on the choice of fc* is 
relatively low. The right plot of figure 2 shows h Tnl (t^ l ) at (S/N) n = 1. Values of r^ L < 10 5 enable to 
detect n TNL > 0.3, again with the PIXIE experiment. 

In the single field case, we can use both the temperatures-distortion correlation C pT or the /x-distortion 
self-correlation C MM to measure n/ NL . As shown in Fig. 3, C pT allows to detect lower values of n/ NL . 



3 Halo bias 



Let us now turn to the effect of running NG parameters onto the halo bias [17,20-22,25,35]. The halo bias 
power spectrum with Gaussian initial conditions can be simply expressed at lowest order in terms of a linear 
(Eulerian) bias parameter 



Ph(k) = (bf) P m (k), 



(29) 



where P m (k) is the dark matter power spectrum. The effect of primordial non-Gaussianity on the halo bias can 
be accurately predicted from a peak-background split [33-37] . As shown in [35] , the non-Gaussian contribution 
to the linear bias induced by a non-zero primordial A-point function is 



Afci(fc) 



A N \Kz) 




din T { s N \k,z 
dine. 



2 Using the explicit expressions above, we find that this condition is verified provided that 



(30) 



b < 10 7 £ 2 . We consider 



the pivots fc* = 0.002 Mpc and fc* = 0.064/iMpc ~ 0.045 Mpc . The former corresponds to the pivot k p of the primordial 
spectrum and the latter to the best pivot value from [26]. For ~ (/nl) 2 ~ anc ^ ^ ~ we nnc ^ that the approximation 
(25) is valid for "/ NL ,n TNL < (0.65 — 0.85) depending on the pivot see Fig. 1 which presents the values of b and b as function of 
the indices n/ NL and n TNlj for the two pivots. One should be aware that in the multiple field case, is larger than ((6/5)/^ L ) 2 , 
so the approximation becomes worse. In general, it seems reasonable to trust our estimation up to ^/ nl i"t nl — 0.5. 



G 




Figure 2: Left: The spectral index n/ NL as function of f^ L at (S/N) nf = 1. Right: The spectral index n rNL as 
function of r^ L at (S/N) nT = 1. Both plots are made for two different pivot scales using fe y = 0.002 Mpc -1 
and n s = 0.96 for PIXIE. ^ 



where 6at are Lagrangian bias parameters, S c ~ 1.68 is the critical threshold for (spherical) collapse and a s is the 
rms variance of the density field at redshift z smoothed on the (small) scale R s of a halo. While this expression 
assumes a universal mass function, it can be generalized to take into account deviations from universality in 
actual halo mass functions [37]. 

The linear matter density contrast 5Az) is related to the curvature perturbation <j>^ during matter domina- 
tion via the Poisson equation. The latter can be expressed as the Fourier space relation SAz) = A4(k, z) 
where 



(31) 



Here, T(k) is the matter transfer function, fi m and H Q are the matter density in critical units and the Hubble 
rate today, and D(z) is the linear growth rate. A4 S is a shorthand for A4(k, z)W(kR s ), where W(kR s ) is a 
spherically symmetric window function (we adopt a top-hat filter throughout this paper). Furthermore, 



1 



4erfP,(k) 



N-2 

n 



d^i 



M s (h,z) 



M s (q, z) 



,k N -2,q,k) 



(32) 



is a projection factor whose /c-dependence is dictated by the exact shape of the A^-point function £^ JV ' ) of 

(3) 

the gravitational potential. For the local constant- /nl model, the factor Ts is equal to /nl in the low k- 
limit (squeezed limit), so that the logarithmic derivative of J-"!^ with respect to the rms variance a s of the 
small-scale density field vanishes on large scales. However, this does not hold for scale-dependent primordial 
non-Gaussianity. In this case, we use expressions (4) and (6) for the bispectrum and trispectrum to evaluate 
the derivative of with respect to a s . 

For generic primordial 3- and 4-point functions, the non-Gaussian halo power spectrum reads 



P h {k) 



E\2 



(6f) 
25 



46f M< 



Tin 



,M) 25 E 
M R {k) 27 1 



hi 



b 2 5 c a 2 R 



1 



din 71 

ciln Or 



riK NL ,M) 
M R (k) 



Pm\k) 



(33) 



9 " x " c M\{k) 

where, on large scales, the last term in the square brackets can generate stochasticity between the halo and 
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Figure 3: Using the temperature- /i-distortion correlation C^ T allows to detect lower values of n/ NL than using 
the //-distortion self-correlation C w in the single field case. The spectral index n/ NL at (S/N) nf is shown for 
two different pivot scales A;*, using k p = 0.002 Mpc" 1 and n s = 0.96 for PIXIE. 



mass density fields if tnl is different from (6/nl/5) 2 [38-41]). We have defined the quantities 

Hnf^M) = \ [ ^q 2 M 2 R (q)P(q)h L (<l), (34) 

Ti(n TOL) M) = 4r / d y o ld L q2 M R (q 1 )M R (q 2 )M R (q 12 )P(q 1 )P(q 2 )r NL (q 1 ,q 12 ), (35) 
°RJ ( 27r ) 

T 2 (n TN „M) = \ [ ^^q 2 1 qlM 2 R (q 1 )M R (q 1 )P(q 1 )P(q^ h (q 1 ,q 2 ) . (36) 

We have used the definitions (8) and (9) to obtain these expressions. We have also emphasized the dependence 
on the parameters i./ NL and « TNL i as well as the halo mass M which, for the top- hat filter, is related to the 
smoothing radius R through R — (3M/47T) 1 / 3 . The values of / NL and at the pivot wavenumber fc* = 0.045 
Mpc -1 are assumed to be known. In the particular case of scale-independent /nl and tnl, i-e. ^/ NL = "-tnl = 0, 
we recover the expressions given in Refs. [38] and [41] . 

In order to assess the ability of forthcoming experiments to probe the scale dependence of the non-linearity 
parameters /nl and tnl through a measurement of the large scale bias, we use the Fisher information content 
on /nl and tnl (see e.g. [17,20-23,25] for application to the scale-dependence of /nl) in the two-point statistics 
of halos and dark matter in Fourier space. 

Computing the Fisher information requires knowledge of the covariance matrix of the halo samples, 

C h (k, M, z) = b 2 (k, M, z)P m {k) + 4, (37) 

n 

where n is the mean number density of the survey. In order to constrain n j NL and n TNL , we assume that we have 
already measured / NL and r NL . Moreover, since we are interested in investigating the possibility of a detection 
of the spectral indices, we take n/ NL = n TNlj = throughout as fiducial values. The Fisher matrix is defined as 
follows 

fdkk 2 1 8C h dC h 

Tij - VWv /sky J ^ 2^1^^' (38) 

where 9i are the parameters whose error we wish to forecast, V smv is the surveyed volume and / s k y is the fraction 
of the sky observed. The integral over the momenta runs from fc m j n = 27r/(V r surv ) 1 / 3 to fc max = 0.03 Mpc _1 //i, 
above which the non-Gaussian bias becomes smaller than contributions from second-order bias and nonlinear 
gravitational evolution. For illustration, we adopt the specifications of a wide-angle, high-redshift survey such 
as BigBOSS or EUCLID: V^ urv / S k y = 50 Gpc 3 //i 3 at median redshift z = 0.7. Furthermore, we ignore redshift 
evolution and assume that all the surveyed volume is at the median redshift. 
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Figure 4: Confidence ellipses obtained by the population of tracers considered with halos with mass larger than 
M = f0 13 M /h , assuming / NL = 20 and 4 = 5x fO 3 (left) and / NL = 50 and 4 = 5x fO 4 (right). 



Table 1 : 1-er errors for the population considered in the two different sets of and r^ L in Fig. 4. 



/nl 


T* 

NL 




NL 


20 


5 x 10 3 


0.30 


0.23 


50 


5 x fO 4 


0.f5 


0.08 



We compute the uncertainties on n / NL and n TNL from a single population of tracers consisting of all halos of 
mass larger than 10 13 M /h. Computing the Lagrangian bias factors from a Sheth-Tormen mass function [42] 
leads a linear and quadratic Lagrangian bias b\ = 0.7 and 62 = —0.4. We take the number density to be 
n = 10- 4 Mpc 3 //i 3 . 

Fig. 4 shows the resulting 68, 95 and 99% confidence contours for the parameters n/ NL and n TNL when we 
assume two different combinations of /j^ L and r^ L . The l-cr errors are displayed in Table 1 . In the specific case 
in which only one degree of freedom is responsible for the perturbations, we can use the relation TNi/fci, kj) = 
^fNh(h)fNh{kj), which leaves us with only one parameter, n/ NL , describing the scale dependence of the 
primordial NG. The 1-er error for n^ NL as a function of is shown in Fig. 5. This result can be compared 
with those of previous work. For a fiducial value of / NL = 50 in particular, we find an error of An/ NL ~ 0.2 
in the case of multi- field models, and Anj NL ~ 0.1 in the case of single-field models. For single-field models, 
this is a factor of 0(3) lower than the forecast error found in Ref. [17] for a survey like EUCLID. We attribute 
this difference to the fact that we have considered the higher-order term 0(/nl) m the halo bias and to the 
parametrization /nl(-?0 = /NL(fc*)(^ 1 ^ 3 /fc*)" /NL considered in Ref. [17] for the running of /nl- In this regards, 
note that K = fciA^&a gives a contribution to the scaling of the external momentum, leading to a suppression 
(for a positive n / NL ) or enhancement (for a negative n y NL ) of the signal with respect to our parametrization in 
Eq. (8). 3 . We have checked that, if we use the parametrization and restrict ourselves to the C(/nl) contribution 
to the halo bias, we are able to reproduce their results. As noted in the introduction, the parametrization used 
in this paper seems to be motivated by various theoretical predictions (see for example [14,25]). 

4 Conclusion 

Even a tiny level of non-Gaussianity in the cosmological perturbations can tell us a lot about the dynamics of 
the inflationary Universe. In this paper, we have focused on local non-Gaussianity, which is a generic prediction 
of multifield inflationary models where cosmological perturbations are sourced by light scalar fields other than 
the inflaton. We have considered the possibility that the non-linear parameter /nl is scale-dependent and, 
extending the previous literature, we have also assumed that tnl may be scale-dependent. This is an unavoidable 
consequence when only a single field other than the inflaton generates the perturbation as the spectral indices 

3 Dctcrmining n ^ NL through /^-distortion using the parametrization of Ref. [17] also leads to a deterioration of the S/N ratio. The 
parameter b is approximated by Eq. (17) with fij- NL replaced by 2n^ NL /3. The correlation C^ T is decreased by a factor of about 
exp(— cn^ NL ) with c ~ 3, 4 for k t = 0.002, 0.045 Mpc -1 respectively, relative to the parametrization Eq.(8). Correspondingly, the 
error a n j is increased by about | exp(cn^ NL ) 
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Figure 5: l-cr error predictions for n/ NL as a function of / N l at the pivot point fc* = 0.045 Mpc -1 for the 
population considered in the case of one-single field models. 

«/ NL and ?v NL are equal. We have considered two possible probes of a running non-Gaussianity. First, we 
have exploited the fact that future measurements of the CMB /i-distortion will be very sensitive to small scales, 
thereby enhancing the effect of a (blue) tilt of the NG parameters. Second, we have assessed the ability of a 
large-scale galaxy survey to constrain the scale dependence of /nl and tnl imprinted in the non-Gaussian halo 
bias. Assuming the detection of a non- vanishing /nl and tnl, we find both for a CMB experiment like PIXIE 
and a large-scale survey like EUCLID that the spectral indices could be measured with an accuracy of O(0.3) 
for /nl = 20 and tnl = 5000. In the case of a measurement of the scale-dependent halo bias, this limit could 
be improved by suitably combining the information from several tracers (e.g. [43]). 
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